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Surface waves, named here as Dyakonov-Tamm waves, can exist at the planar interface of 
an isotropic dielectric material and a chiral sculptured thin film (STF). Due to the periodic 
nonhomogeneity of a chiral STF, the range of the refractive index of the isotropic material is 
smaller but the range of the propagation direction in the interface plane is much larger, in 
comparison to those for the existence of Dyakonov waves at the planar interface of an isotropic 
dielectric material and a columnar thin film. 
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1 INTRODUCTION 

Less than two decades ago, Dyakonov [1] theoretically predicted the propagation of a surface 
wave at the planar interface of an isotropic dielectric material and a positively uniaxial dielectric 
material with its optic axis wholly parallel to the interface plane. If ip indicates the angle between 
the optic axis and the direction of surface-wave propagation, and n s is the refractive index of 
the isotropic dielectric material, then the Dyakonov wave exists for rather narrow ranges of ip 
and n s . The consequent significance of Dyakonov waves for optical sensing and waveguiding was 
recognized thereafter [2,3]. Since then, the concept of the Dyakonov wave has been extended 
to the planar interfaces of isotropic and biaxial dielectric materials [4]. The possibility of the 
anisotropic material being artificially engineered, either as a photonic crystal with a short period 
in comparison to the wavelength [5] or as a columnar thin film (CTF) [6], has also emerged. Let 
us note here that the Dyakonov wave still remains to be experimentally observed, in part due 
to the narrow range of ip for its existence [5]. 

The anisotropic material is taken to be homogeneous in all of the foregoing and other reports on 
the Dyakonov wave. What if the anisotropic material were to be chosen as periodically nonho- 
mogeneous in a direction normal to the bimaterial interface? This question initiated a research 




Figure 1: Geometry of a structurally right-handed helix. 



project, the first results of which are being communicated here. Being a natural extension of 
a CTF, a chiral sculptured thin film (STF) was chosen as the periodically nonhomogeneous 
anisotropic material [7]. 

A chiral STF is made by directing a vapor flux in vacuum at an oblique angle onto a rotating 
substrate. Under suitable conditions, an assembly of parallel nanohelixes of the evaporated 
species forms, with the helical axes perpendicular to the substrate. An example of a single 
nanohelix is illustrated in Figure [TJ By adjusting deposition parameters, both the pitch 2Q 
and the angle of inclination \ G (0, vr/2] can be controlled. Each nanohelix, composed of 
multimolecular clusters with ~ 3 nm diameter, is effectively a continuously bent column of 
~ 100-nm cross-sectional diameter. Therefore, at visible frequencies and lower, a chiral STF 
may be regarded as a linear, locally orthorhombic, unidirectionally nonhomogeneous continuum 
whose relative permittivity dyadic is akin to that of chiral smectic liquid crystals [8] . 

In formulating the surface-wave-propagation problem on the planar interface of an isotropic, 
homogeneous, dielectric material and a chiral STF, we adopted a methodology originally devel- 
oped by Tamm in 1932 for a realistic Kronig-Penney model. Instead of assuming the solid to 
occupy the entire space, as is commonplace in solid-state physics [9], Tamm assumed the solid to 
occupy only a half-space. The incorporation of the oft-neglected surface led to the emergence of 
electronic states localized to the surface. Tamm states were experimentally observed in 1990 on 
the surfaces of superlattices [10], and their optical analogs for superlattices of isotropic materials 
are being investigated these days [11,12]. 

Given the braiding of Dyakonov waves and Tamm states in this communication, we decided to 
name the surface wave at the planar interface of an isotropic, homogeneous, dielectric material 
and a chiral STF as the Dyakonov-Tamm wave. Section [2] presents the boundary-value prob- 
lem and the dispersion equation for the Dyakonov-Tamm wave. Section [3] contains numerical 
results when the chiral STF is chosen to be made of titanium oxide [6, 13]. An exp(— iuit) time- 
dependence is implicit, with u denoting the angular frequency. The free-space wavenumber, the 
free-space wavelength, and the intrinsic impedance of free space are denoted by k Q = tj^/e /i , 
A = 2-K/k , and r\ Q = \J fJ, /^o, respectively, with {i Q and e D being the permeability and per- 
mittivity of free space. Vectors are underlined, dyadics underlined twice; column vectors are 
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underlined and enclosed within square brackets, while matrixes are underlined twice and simi- 
larly bracketed. Cartesian unit vectors are identified as u x , u y and u z . 

2 FORMULATION 

2.1 Geometry and permittivity 

Let the half-space z < be occupied by an isotropic, homogeneous, nondissipative, dielectric 
material of refractive index n s . The region z > is occupied by a chiral STF with unidirectionally 
nonhomogeneous permittivity dyadic given by [7] 

g(z) = e g z (z)-S y (x)-g re{ -gy(x)-g(z), z>0, (1) 

where the reference relative permittivity dyadic 

g ref = e a u z u z + e b u x u x + e c u y u y (2) 

indicates the locally orthorhombic symmetry of the chiral STF. The dyadic function 

gjjS) = COS + (U^UJ, + UyUy) 

(HZ \ 
+ tp) [UyUv - u x u y ) + u z u z (3) 

contains 20 as the structural period, ip as an angular offset, and h = ±1 as the handedness 
parameter. The tilt dyadic 

Sy(x) = (MA + u z u z ) cos X 

+ {u z u x - u^u z ) sin x + u y u y (4) 

involves the angle of inclination x- The superscript T denotes the transpose. 

Without loss of generality, we take the Dyakonov-Tamm wave to propagate parallel to the x axis 
in the plane z = 0. There is no dependence on the y coordinate, whereas the Dyakonov-Tamm 
wave must attenuate as z — > ±oo. 

2.2 Field representations 

In the region z < 0, the wave vector may be written as 

k s = ku x — a s u z , (5) 

where 

k 2 + a 2 s = k 2 n 2 s , (6) 

k is positive and real-valued for unattenuated propagation along the x axis, and Im [a s ] > for 
attenuation as z — ► — oo. Accordingly, the field phasors in the region z < may be written as 



E(r) 



A 1 u y + A 2 [—u x + —u, 



exp(ik s -r), z<0, (7) 
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and 



K(r) = r,- 1 



A\ — u„ + — u. 



k k n 
where A\ and A 2 are unknown scalars. 



A 2 n 2 s u y 



exp(ik s • r ) , z < , 



The field representation in the region z > is more complicated. It is best to write 

E(r) = e(z) exp(ira) 1 

> , z > . 

H(r) = h(z) ex.p(inx) 

and create the column vector 

[/(*)] = \ e x{z) e y {z) h x {z) h y (z)] T . 
This column vector satisfies the matrix differential equation [7] 



fl r 7T 7 1 

^ [/(*)] = i\E& + f/>, k) I • [m] , z>o. 



where the 4x4 matrix 
£(C,«)] 






/i e D (e c - ed) cos C sin ( 
e Q (e c sin 2 ( + e d cos 2 C) 



/i D 

-fJL 



-e (e c cos 2 C + £d sin 2 C) 
-he Q (e c — erf) cos£ sin£ 



+ K sm x cos x 



--J51-&L 



cos C h sin £ 












—h sin£ 
cos C 



+ 


















and 



Q/ 



e a cos 2 x + e b sin 2 x 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



Two independent techniques [14,15] exist to solve (|lip . which may be harnessed to determine 
the matrix [N] that appears in the relation 



mm = [N\ . [/( +)] 



(14) 



to characterize the optical response of one period of the chiral STF. By virtue of the Floquet- 
Lyapunov theorem [16], we can define a matrix [Q] such that 



[K] =exp|i20[Q]} 



(15) 
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Both [N] and [Q] share the same eigenvectors, and their eigenvalues are also related. Let [t}^ n \ 
(n = 1,2,3,4), be the eigenvector corresponding to the nth eigenvalue a n of [N] \ then, the 
corresponding eigenvalue a n of [Q] is given by 



a. 



.In a n 



(16) 



2.3 Dispersion equation for Dyakonov Tamm wave 

For the Dyakonov-Tamm wave to propagate along the x axis, we must ensure that Imfai^] > 0, 
and set 



[/(0+)]= [t] (2) 



Bo 



(17) 



where B\ and B2 are unknown scalars; the other two eigenvalues of [Q] describe waves that 
amplify as z — > 00 and cannot therefore contribute to the Dyakonov-Tamm wave. At the same 
time, 



[/(o- 








1 





as „-l 








2 — 1 

-n% r) 



(18) 



by virtue of ([7]) and ((8j). Continuity of the tangential components of the electric and magnetic 
field phasors across the plane z = requires that 



[/(o-)] = [/(0+)] 



(19) 



which may be rearranged as 



" At ' 







A 2 







Bi 







B 2 








[M] 



For a nontrivial solution, the 4x4 matrix [M] must be singular, so that 

det [K] = 

is the dispersion equation for the Dyakonov-Tamm wave. 



(20) 



(21) 



3 NUMERICAL RESULTS AND DISCUSSION 

Although chiral STFs may be made by evaporating a wide variety of materials [7, Chap. 1], 
the constitutive parameters of chiral STFs have not been extensively measured. However, the 
constitutive parameters of certain columnar thin films (CTFs) are known. CTFs are assemblies 
of nanorods oriented at an angle % to the substrate and are produced by directing the vapor 
at an angle Xv onto a stationary substrate, as shown in Fig. [21 the vapor incidence angle 
Xv ( m addition to the evaporant species) determines the constitutive parameters [13]. When 
the substrate is rotated about a normal passing through it at a constant angular velocity of 
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Figure 2: Schematic of the growth of a columnar thin film. The vapor flux is directed at an 
angle Xv, whereas nanorods grow at an angle x > Xv 



reasonable magnitude, parallel nanohelixes grow instead of parallel nanorods, and a chiral STF 
is deposited instead of a CTF [7,17]. Although the substrate is nonstationary, the functional 
relationships connecting e a ^, c and x to Xv for CTFs would substantially apply for chiral STFs, 
since the vapor incidence angle Xv remains constant during the deposition of thin films of either 
kind. 

Among the CTFs which have been characterized are those made of titanium oxide, a material 
important in many practical applications [6]. Empirical relationships have been determined for 
titanium-oxide CTFs at A D = 633 nm by Hodgkinson et al. [13] as 



1.0443 + 2.7394 



1.6765 + 1.5649 



1.3586 + 2.1109 



Xv 

vr/2 

Xv 

tt/2 

Xv 

tt/2 



1.3697 



0.7825 



1.0554 



Xv 

tt/2 

Xv 

tt/2 

Xv 

tt/2 



and 



tanx 



2.8818 tan X v , 



(22) 
(23) 
(24) 

(25) 



where Xv arid x are in radians. We must caution that the foregoing expressions are applicable 
to CTFs produced by one particular experimental apparatus, but may have to be modified for 
CTFs produced by others on different apparatuses; hence, we used these expressions for the 
numerical results presented in this section for chiral STFs simply for illustration. Furthermore, 
we set h = 1, O = 197 nm, and Xv — 7.2°. Following Walker et al. [4] and Polo et al. [6], we left 
tp and n s as variable parameters. All numerical results presented in this section were computed 
for A = 633 nm. 

As mentioned in Section^ the matrix [N] can be calculated using two numerical techniques [15]: 
(i) the piecewise uniform approximation technique and a series technique based on the Maclaurin 
expansion of [P(C> k)] with respect to £. Both yield the same results, and the piecewise uniform 
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Figure 3: v as a function of ip with Xv = 7.2° for n a = 1.631, 1.635, 1.64, 1.645, and 1.65. 

approximation technique was selected for calculations reported here. Basically, the technique 
consists of subdividing the chiral STF into a series of electrically thin sublayers parallel to the 
interface, and assuming the dielectric properties to be spatially uniform in each sublayer. The 
accuracy of this technique depends on the thickness of the sublayers, with thinner ones yielding 
more accurate results. Based on experience [7,15], a sublayer thickness of 2 nm gives reasonable 
results. 

The magnitude of the phase velocity of the Dyakonov-Tamm wave was compared with that of 
the phase velocity of the electromagnetic wave in the bulk isotropic material. For this purpose, 
we defined the relative phase speed 

v = v DT /v s , (26) 

where vdt = w/k is the phase speed of the Dyakonov-Tamm wave and v s = \ jn sy Je [i is the 
phase speed of the electromagnetic wave in the bulk isotropic material. Figure [3] shows v as a 
function of tj; for several values of n s . The phase velocity of the Dyakonov-Tamm wave, like 
several other surface waves [1]- [6], was found to be lower in magnitude than the phase velocity 
of the electromagnetic wave in the bulk isotropic material. 

The minimum and maximum values of n s (1.631 and 1.65, respectively) in Figure [3] represent 
the approximate limits of the ra s -range for which the determinantal equation (|21D representing 
the boundary conditions between the two material could be solved. Outside this n s -range, the 
Dyakonov-Tamm wave can not exist for the chosen parameters. 

Each of the curves in Figure [3] was drawn over the continuous ?/;-range for which the Dyakonov- 
Tamm wave was found to exist. The plot in Figure [3]is restricted to tp G [—20°, 100°]. For each 
value of tp in this plot (and all remaining plots in this paper), a similar point exists at ip + 180° 
with an identical value of the function; thus, for each value of n s , there exist two separate ranges 
of ip G [—180°, 180°] over which the Dyakonov-Tamm waves exist. Of the curves presented, the 
one for n s = 1.4, about mid-range in n s , has the widest range in ip with a width Aip = 98°. 
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The mid-point of the ^-range is ip m = 37°. As n s approaches either end of the n^-range, 
diminishes. At the low end of the n s -range, Atp = 51° when n s = 1.635, but Atp = 20° when 
n s = 1.631. At the high end of the n s range, Atp = 54° when n s = 1.645, but Aip = 15° 
when n s = 1.65. Only a slight variation in tp m , the mid-point of the ip-r&nge, is seen. For 
n s = 1.631, 1.635, 1.4, 1.645, and 1.65, we find ip m = 36°, 36.5°, 37°, 40°, and 39.5°, respectively. 
There seems to be a slight increase in ip m as n s increases. The values of ip m , however, are only 
approximate since the end-points of the ■i/'-range were only determined by the last whole degree 
lying inside the range. This may account for why ip m at n s = 1.65 is slightly lower than that at 
n s = 1.645. 

Every curve in Figure is smooth with a broad minimum in the vicinity of 35° to 40° which 
levels off at both ends of the ^-range. The minimum is deepest for curves representing mid- 
range values of n s , while curves at extreme values of n s are nearly flat. As n s decreases, the v 
vs. ip curve shifts downward. 

The confinement of the Dyakonov-Tamm wave to the interface is described by the decay con- 
stants which are given by the imaginary part of the two eigenvalues in the chiral STF (a± and o^) 
and the single eigenvalue a s in the isotropic dielectric material. We found all three eigenvalues 
to be purely imaginary, which are shown in Figure [Has functions of ip for several values of n s . 

In Figure [U typical values of Im[«i] are about one order of magnitude larger than Imf^]. Every 
Im[ai] vs. ip curve in Figure Hh, is bell-shaped with a maximum in the vicinity of ip = 35° to 
40°. The maximum change in Im[ai] over the ■(/'-range is greatest for n s = 1.64, the mid-range 
of n s . At this value of n s , the variation in Im[«i] is still less than 0.5%. Toward each end of the 
n s -range, the Im[ai] vs. ip curve flattens. 

The Im[a2] vs. ip curves in Figure0b show two types of behavior. At the high end of the ra s -range 
(n s = 1.645,1.65), the curves are bell-shaped with greater variation occurring at n s = 1.645 
away from the end of the ra^-range, just as evinced by Im[ai]. In contrast, in the lower half of 
the ra s -range (n s = 1.631, 1.635, 1.64), the Imfo^] vs. ip curves do not level off at the ends of the 
respective ^-ranges, but continue to decrease all the way to zero. 

The behavior of the decay constant in the isotropic material can be gleaned from the Im[a s ] 
vs. ip curves in Figure Eh. These curves look very similar to those of Imfc^] in Figure Hb. 
However, the dependence of the shape of the Im[a s ] vs. ip curve on n s is opposite that of the 
Im[o!2] vs - ip curve. The curves in Figure 0b for low values of n s are bell-shaped with Im[a s ] 
remaining non-zero over the entire ^-range; but the curves for large values of n s maintain their 
downward curvature with Im[a s ] going to zero at both ends of the respective ^-ranges. Thus, 
the Dyakonov-Tamm wave becomes delocalized from the interface z = at the limiting values 
of n s ; derealization occurs on the chiral STF side of the interface at low values of n s , but on 
the isotropic substrate side of the interface for high values of n s . 

Similar results are obtained at other values of Xv As an example, Figure [5] displays v when 
Xv = 25°. As in Figure El the maximum and minimum values of n s displayed in Figure [5] 
represent the approximate limits of the range of n s over which surface-wave propagation is 
possible. The width of the n s -range for Xv = 25° is roughly half of that obtained for Xv = 7.2°. 
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Figure 4: Decay constants as functions of ip at A D = 633 nm for same values of n s and \v as in 
Figure [3j a) a\ b) a<i c) a s . 
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Figure 5: v as a function of ip with \v = 25° for n s = 1.947, 1.95, 1.955, 1.958. 

The curves for Xv = 25° display the same general shape and trends as for the case of Xv = 7.2°, 
with similar values of Atp. However, the minimums have shifted to a lower value of tp between 
15° and 20°. Now, tp m = 13°, 14.5°, 17.5°, and 18° for n s = 1.947, 1.950, 1.955, and 1.958, 
respectively. Again, a slight upward drift of ip m as n s increases can be seen. The figure also 
shows that the minimum value of v is larger for Xv = 25° than it is for Xv = 7.2°. Although not 
shown, all three decay constants are lower at the higher value of Xv- 

Calculations were also made at A G = 533 and 733 nm, with the assumption of the same consti- 
tutive and geometric parameters as at A D = 633 nm. Qualitatively similar conclusions on the 
width of the ranges of n s and ijj for Dyakonov-Tamm waves were drawn from the numerical 
results obtained. Furthermore, the minimum value of v was found to increase and the maximum 
values of Im[ai], Im[o2], and Im[ct s ] were found to decrease, as the ratio A /fi increases. 

Additional calculations, not shown here, were made for an identical chiral STF except with 
structural left-handedness {h = —1). The exact same results were obtained with the left- 
handed chiral STF as with the right-handed chiral STF presented in Figures 2-4. It must borne 
in mind that, by changing h from ±1 to =Fl in we invert not only the structural handedness 
of the chiral STF but also the sense of rotation brought about by ijj ^ 0. 



To conclude, we examined the phenomenon of surface-wave propagation at the planar interface 
of an isotropic dielectric material and a chiral sculptured thin film. The boundary- value problem 
was formulated by marrying the usual formalism for the Dyakonov wave at the planar interface 
of an isotropic dielectric material and a columnar thin film with the methodology for Tamm 
states in solid-state physics. The solution of the boundary-value problem let us deduce the 
existence of Dyakonov-Tamm waves. 

In constitutive terms, the major difference between a CTF and a chiral STF is the periodic 



4 CONCLUDING REMARKS 
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nonhomogeneity enshrined in S. (z); in the limit — > oo, a nanohelix uncurls into a nanorod, 
and a chiral STF transmutes into a CTF. However, the distinction between chiral STFs and CTFs 
is nontrivial, as may be deduced from the Floquet-Lyapunov theorem [16, 18]. In comparison 
to the Dyakonov wave localized to the planar interface of an isotropic dielectric material and a 
CTF [6] , we found that the n s -range for the existence of a Dyakonov-Tamm wave at the planar 
interface of an isotropic dielectric material and a chiral STF is smaller. However, the ip-vange is 
much larger in width: in comparison to Atp < 1° with CTFs [6] — and < 5° with effectively 
uniaxial, short-period photonic crystals [5] - the width of the ^-range is as high as 98° in 
Figure [3] with chiral STFs. This implies that Dyakonov-Tamm waves could be detected much 
more easily than Dyakonov waves. 

Acknowledgment. A.L. gratefully acknowledges a fascinating discussion with Profs. J. Mar- 
torell (Universidad de Barcelona) and D.W.L. Sprung (McMaster University). 
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